Torqued fireballs in relativistic heavy-ion collisions 
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We show that the fluctuations in the wounded-nucleon model of the initial stage of relativistic 
heavy-ion collisions, together with the natural assumption that the forward (backward) moving 
wounded nucleons emit particles preferably in the forward (backward) direction, lead to an event- 
by-event torqued fireball. The principal axes associated with the transverse shape are rotated 
in the forward region in the opposite direction than in the backward region. On the average, 
the standard deviation of the relative torque angle between the forward and backward rapidity 
regions is about 20° for the central and 10° for the mid-peripheral collisions. The hydrodynamic 
expansion of a torqued fireball leads to a torqued collective flow, yielding, in turn, torqued principal 
axes of the transverse-momentum distributions at different rapidities. We propose experimental 
measures, based on cumulants involving particles in different rapidity regions, which should allow 
for a quantitative determination of the effect from the data. To estimate the non-flow contributions 
from resonance decays we run Monte Carlo simulations with THERMINATOR. If the event-by-event 
torque effect is found in the data, it will support the assumptions concerning the fluctuations in the 
early stage of the fireball formation, as well as the hypothesis of the asymmetric rapidity shape of 
the emission functions of the moving sources in the nucleus-nucleus collisions. 

PACS numbers: 25.75.-q, 25.75Gz, 25.75.Ld 
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I. INTRODUCTION 

It is believed that the forward-backward (FB) rapid- 
ity correlations may reveal important information on 
the mechanism of particle production in high-energy 
hadronic and nuclear collisions. Long-range rapidity cor- 
relations uncover properties of the dynamical system at 
a very early stage. For that reason the FB multiplic- 
ity correlations have been studied experimentally 0, 
and theoretically In this paper we show that the 

wounded-nucleon approach [9j to the initial stage of the 
collisions leads to a new manifestation of the FB fluc- 
tuations: the torqued fireball. Specifically, parts of the 
fireball are rotated in the transverse plane in one direc- 
tion in the forward rapidity region, and in the opposite 
direction in the backward rapidity region. The magni- 
tude and sign of the torque angle fluctuate from event to 
event. 

The following ingredients are responsible for the ap- 
pearance of the torque effect: 

1. statistical fluctuations of the transverse density of 
the sources (wounded nucleons) [To|, and 

2. the asymmetric shape of the particle emission 
function, peaked in the forward (backward) rapid- 
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ity for the forward (backward) moving wounded nu- 
cleons. 

On the average, for the studied Au+Au collisions at 
the highest energies available at the BNL Relativistic 
Heavy Ion Collider (RHIC), the relative torque angle be- 
tween the principal axes in the forward and backward 
rapidity regions is about 20° for the central, and 10° for 
the mid-peripheral collisions. On the other hand, an ex- 
perimental observation of nonzero torque angles of the 
expanding fireball could shed light on the mechanism of 
the formation of dense matter. A finite torque between 
reaction planes at different rapidities could influence el- 
liptic and directed flow studies using reaction planes de- 
termined in different pseudorapidity intervals [11| . 

The paper has three parts, referring subsequently 
to the early wounded-nucleon phase of the colli- 
sion (Sect. [n|, the intermediate hydrodynamic stage 
(Sect. IIII[) . and, finally, the statistical hadronization 
phase (Sect. IIV[) . where hadrons are produced. In 
Sect. Ill Al we introduce the necessary elements of the 
wounded-nucleon approach, in particular, the asymmet- 
ric rapidity-dependent emission functions 04^1 of the 
forward and backward moving wounded nucleons. In 
Sect. Ill Bl we explain how statistical fluctuations of the 
density of the forward and backward moving wounded 
nucleons in the transverse plane lead to the event-by- 
event torque effect. Simulations are carried out with 
GLISSANDO in the so-called mixed model 0, [iljj . 
incorporating an admixture of binary collisions into the 
wounded-nucleon model. We introduce quantitative 
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measures of the torque in Sec. Ill CI Specifically, we use 
the difference of the forward and backward angles of the 
principal axes, as well as three angles, forward, back- 
ward, and central, to construct useful statistical quanti- 
ties. Next, in Sect. IIIII we present the results of running 
the 3+1 dimensional (perfect fluid) hydrodynamics, as 
described in [l5j]. The calculations show that the torque 
survives the hydrodynamic stage of the evolution. Then, 
in Sect. HV Al we pass to discussing the final stage, namely, 
the statistical hadronization (for a review, see, e.g. [la]). 
This stage would wash out the torque effect, unless care- 
ful experimental measures are used. This is because even 
at a fixed spatial geometry of the fireball the finiteness 
of the number of produced particles causes large fluctu- 
ations of the event-plane angle, covering up the torque 
angle. We thus propose to investigate measures based on 
cumulants [l7j . introduced in Sect. IIV Bl They are con- 
structed in such a way that the relative FB torque angle 
can be extracted. We examine the non-flow contribu- 
tions to these measures via simulations in THERMINATOR 
[l8j . proving that it is possible to find the torque effect 
with the proposed methods in the large-statistics RHIC 
data. In Sect. IIVDI we repeat this analysis for the par- 
ticles taken from three rapidity bins, forward, backward, 
and central, with similar conclusions. 



II. TORQUED FIREBALL 

In this Section we describe the earliest stage of the 
relativistic heavy-ion collision in terms of the wounded- 
nucleon model |!) . All simulations are carried out for the 
Au+Au collisions at the highest RHIC energy of ^snn = 
200 GeV. 



A. Wounded nucleons with rapidity profiles 

The wounded-nucleon Q approach is commonly used 
to describe the early stage of the heavy-ion collisions. In 
the collision process N w nucleons get wounded, i.e. inter- 
act inelastically at least once, as well as JVbin binary colli- 
sions occur. Roughly speaking, wounded nucleon (partic- 
ipants) are responsible for the soft emission, while binary 
collisions describe hard processes. Both serve as sources 
for the formation of the density of energy or entropy in 
the initial fireball. At RHIC, the mixed model [l3|, where 
the total number of the produced particles is given by the 
combination 

^prod = A fl^-N w + aN hi ^\ , (1) 

describes quite successfully the multiplicity data flij . 
The proportionality constant A depends on the energy 
of the collision but not on its centrality. Throughout this 
work we use a = 0.145, corresponding to collisions at 
y/sNN = 200 GeV [lH- 




FIG. 1. (Color online) The emission profiles in space-time ra- 
pidity for the wounded nucleons (dashed lines) and the binary 
collisions (solid line). The profile /+ (/_) corresponds to the 
forward (backward) moving wounded nucleons. 



The concept of sources can be extended to account for 
the rapidity dependence of the produced particles. The 
pseudorapidity distribution is given as a sum of contri- 
butions from the emission of forward and backward go- 
ing wounded nucleons. Within such a scheme Bialas and 
Czyz successfully described [4] the distribution of charged 
particles in pseudorapidity in the deuteron-gold collisions 
lSj. The independent emission from the forward and 
backward going nucleons determines specific FB multi- 
plicity correlations [1, 0] ■ in particular, this hypothesis 
has been tested successfully in Ref. Q with the FB mul- 
tiplicity correlation data from the PHOBOS collabora- 
tion 0|. Based on this idea, Ref. [§[ postulated that 
in nucleus-nucleus collisions the emission profile defining 
the initial density (in the space-time rapidity rj\\ and the 
transverse plane coordinates x, y) has the form 

F(v\\,x,y) = 0--a)[p+{x,y)f+(r]\\)+P-(x,y)f-{r]i\)] 
+ apbin{x,y)f(ri\\), (2) 

where p±(x,y) is the density of the forward and back- 
ward going wounded nucleons at a given point in the 
transverse plane, Pbm{x,y) is the binary collisions den- 
sity, while /+(f7||) and f-{r]\\) describe the corresponding 
wounded-nucleon emission profiles. Finally, f(j)\\) is the 
emission profile for the binary collisions. These func- 
tions are chosen appropriately, taking into account the 
following features: the profile /+ (/_) is peaked in the 
forward (backward) direction, i.e. the wounded nucleon 
emits mostly in its own forward hemisphere, with a quite 
broad distribution, whereas the binary profile / is (for 
the collision of identical nuclei) symmetric. 

We remark that the asymmetric wounded-nucleon 
emission functions in Eq. ^ give a tilt away from the 
beam axis for the initial fireball (in the (x,r]\\) plane). 
The hydrodynamic expansion Q of such a tilted fireball 
generates the directed flow and for the first time the ex- 
perimental observations at ^/snn = 200 GeV 20] could 
be reproduced. 
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FIG. 2. A cartoon visualization of the torque effect. A ran- 
dom cluster of wounded nucleons, drawn here at the edge of 
the ellipse, with three nucleons moving in the forward (F) 
direction (open circles) one one moving in the backward (B) 
direction (filled circle), causes a random torque of the prin- 
cipal axes. The angle of the torque is higher in the forward 
direction than in the backward direction. In the the central 
rapidity region (C) the effect is between the F and B cases. 



Following Ref. @|. 
izations: 

f(v\\) = exp 



we choose the following parameter- 
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-0{\V\\\ ~ Vo. 



Mv\\) = fF(v\\)f(v\\), 
f-(V\\) = fF(-V\\)f(V\\), 



with 



fF(V\\) 



2%, 



0, 



1, 



V\\ < -Vm 
-Vm < V\\ < Vrn 

Vrn < V\\ 



(3) 



(4) 



The values of parameters, describing the RHIC data after 
the hydrodynamic evolution [8|, are 



m = i) 

r\ m = 3.36, 
a n = 1.3. 



(5) 



The profile functions are shown in Fig.[TJ We note that by 
construction / + (?7||) + f-{rj^) = /(?7||). Parametrization 
(U]) is chosen in such a way, that after the hydrodynamic 
evolution and statistical emission 0, [l5[ one correctly 
reproduces the spectra [2l[ of particles produced at dif- 
ferent rapidities in the Au+Au collisions at the highest 
RHIC energy. 



B. Fluctuations and the torque 

The initial density of the fireball ([2]) can be obtained 
in a Glauber Monte Carlo approach. The densities of the 
forward and backward going wounded nucleons p±(x,y) 
and of the binary collisions Pbin(x,y) are obtained with 



FIG. 3. (Color online) The schematic figure of the torqued 
fireball, elongated along the rj« axis. The direction of the 
principal axes in the transverse plane rotates as rj^ increases. 
The left and right pictures correspond to the rank-2 (elliptic) 
and rank-3 (triangular) cases, respectively. The effect occurs 
event-by-event . 



Monte Carlo simulations by GLISSANDO 12J. These dis- 
tributions fluctuate on event-by-event basis. The phe- 
nomenon has a purely statistical origin, as the positions 
on nucleons in nuclei fluctuate. The event- by- event fluc- 
tuations in the wounded-nucleon approach are know to 
cause important effects, such as the increase of the ellip- 
tic deformation, resulting in larger elliptic flow [loj . or 
the recently discussed triangular flow [221 ] - 

Another effect due to fluctuations, focal to this work, 
is the event-by-event torque of the fireball. Its appear- 
ance is simple to understand. For the sake of simplicity 
let us consider the situation depicted in Fig. O A clus- 
ter of wounded nucleons is formed, here drawn at the 
edge of the fireball. It contains three wounded nucleons 
moving forward and one moving backward. The cluster 
causes the twist of the principal axis. However, due to 
the shape of the emission functions of Fig. Q] the shift is 
different at various values of the space-time rapidity. At 
forward t]\\ ~ 3 there are practically no backward moving 
wounded nucleons, hence the three forward-moving nu- 
cleons cause the torque. In the backward direction only 
one wounded nucleon from the cluster causes the torque, 
while in the central r)\\ region all four nucleons contribute, 
but according to Eq. (0| their relative weight is reduced 
be a factor of 2. Thus the relative weight of the effect of 
the cluster in the forward, backward, and central regions 
is as 3 : 1 : 2. The result is the schematic arrangement of 
the principal axes as drawn in Fig. [5] 

In an actual Monte Carlo simulation many clusters oc- 
cur and the situation is more complicated, but the origin 
is as described above. The effect appears on the event- 
by-event basis and by symmetry the mean value of the 
torque angle vanishes upon averaging over events. Thus, 
the effect may only be revealed in event-by-event studies 
of fluctuations, see Sect. Ill CI 

A similar phenomenon occurs for the axes of the trian- 
gular shape and the axes corresponding to higher Fourier 
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moments in the azimuthal angle. The shape of the fire- 
ball is depicted schematically in Fig. [31 where the torque 
angle, somewhat exaggerated, is shown for the elliptic 
and triangular deformations. As mentioned above, the 
torque appears on the event-by-event basis, varying in 
the direction and value. 



C. Characteristics of the torqued the fireball 

In this subsection we provide quantitative studies of 
the torque effect at the instant of formation. In a given 
event, the angle of the principal axis for the Fourier mo- 
ment of rank-fc for n sources at some 7711 is given by the 
formula 



^( fc ) 



I arctan ( Eti sm(k^) 
k \J2i=l W i r i cos(fo^) 



(6) 



where (fj, 4>i) are the polar coordinates of the source po- 
sition with respect to the center of mass of the slice of 
the fireball at some fixed 7711, and Wi is the weight of 
the source. Explicitly, for the forward-moving wounded 
nucleons w; = (1 — a)/+(?7||), for the backward ones 
Wi = (1 — a)/_(?7||), while for the binary collisions 
All spatial angles are measured relative to 
the axes perpendicular to the reaction plane. The inter- 
pretation of the angle \£W is that the azimuthal Fourier 
moment of rank k is highest along that direction. The 
angle is defined modulo 2n/k. For k — 2 the angle ty( 2 ' is 
the angle of the principal axes of the moment of inertia. 
For brevity of notation, we shall skip the superscript (2) 
from the rank-2 quantities, while retaining superscripts 
indicating higher ranks. 

The simplest measure of the torque effect is the differ- 
ence of the angles at forward and backward values 

of 7711, 



AWfa||) = *< fc >fa||)-* 



(fc) 



(7) 



As argued above, this quantity fluctuates event-by-event. 
The result of the GLISSANDO simulations is shown in 
Fig. |U We plot the event-by-event distribution of Aps 
for the 20-30% centrality class and 77|| = 0.5 (the narrower 
distribution) and m — 2.5 (the wider distribution). Ac- 
cording to the previous qualitative discussion illustrated 
with Fig. [21 the widening of the distribution with increas- 
ing 7711 is expected. The rms width of the distributions 
of Fig. H for c = 20 - 30% are 2° for = 0.5 and 7° 
for 7711 = 2.5. An analogous study for c = — 10° (not 
shown) yields 4° for 7711 = 0.5 and 16° for 7711 = 2.5. Thus 
the spread enhances with increasing 7711 , moreover, the 
distributions are wider for the central collisions. 

For the rank-3 analysis (triangular flow) we have qual- 
itatively the same effect. For c = 20 — 30% the rms width 
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A FB [deg] 

FIG. 4. Distribution of the difference of the forward and 
backward torque angles, ^f — ^b, for the elliptic deformation. 
The narrower and wider distributions correspond to the space- 
time rapidities 7711 = 0.5 and 2.5, respectively. Centrality 20 — 
30%, mixed model for Au+Au collisions, a — 0.145. 
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One could also overlay a statistical distribution of weights, as 
was done in Ref. [13 . 



FIG. 5. The dependence of the rms width of the Afb distri- 
bution on ?7|| at various centralities c. 



are 4° for 7711 = 0.5 and 16° for 7711 = 2.5, which are larger 
than for the rank-2 case at the same centrality. 

In Figs. [5] and [5] we show the dependence of the rms 
width of the distributions of and Ap B as functions 
of 7711 at various centralities. For the rank-2 case we note 
that the widths grow gradually from at 7711 = to 10 — 
20° at 7711 = 4, with the largest angle for the most central 
events. The rank-3 widths are somewhat larger than in 
the rank-2 case, except for the most central case. We 
note that at fixed 7711 the dependence on centrality is non- 
monotonic, with the lowest value for c = 20 — 30% for the 
rank-2 case and 40 — 50% for the rank-3 case. We also 
note that for the rank-2 case the most central collisions 
lead to significantly larger torque fluctuations than for 
the less central collisions, cf. Fig. [SI 
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FIG. 6. Same as Fig. [5] for the rank-3 case (triangular flow). 

One may introduce other measures, involving three an- 
gles associated with the forward, backward, and central 
77|| regions. This might be advantageous, as the central 
rapidity region is experimentally better covered experi- 
mentally. We define the two relative angles 

=^( fe )(?7||)-*( fe )(0), 
A<& = *W(-»7||)-*< fc )(0), (8) 
and their covariance, 

cov fbc = (^FC^Bc) evcnts - (^) 

The forward-backward-central correlation coefficient is 
defined as 

P P BC ee cov^ c /(«r(A« MA<& )). (10) 

In Fig. [7] we present the 2-dimensional distribution plot 
of Ape and A^c for the rank-2 case for m = ±2.5 and 
the 50-60% centrality class-where the anti-correlation of 
the angles is clearly visibleo 

Figures [5] and |H] show the rapidity dependence of 
cov fbc ^ or * ne ran k-2 and rank-3 cases. We note that 
these measures drop monotonically with 7711 from zero to 
negative values in the range —0.005 to —0.02, reaching a 
plateau near 7711 ~ 3.5. For a fixed r]\\ the dependence on 
the centrality class is non-monotonic. From the data of 
Figs. [5] and [5] one may obtain the correlation of Eq. (ITU1) . 
This quantity grows from the value — 1 at 7711 =0 up to 
about —0.5 at r]\\ = 4, similarly for the rank-2 and rank-3 
cases. 



2 We note that the points in Fig. 0occupy all quadrants, i.e. there 
are cases where both shifts Ape and Age have the same sign. 
This is because for each rapidity we evaluate independently the 
center of mass, which is the origin for the transverse coordinate 
system. If we evaluated the angle in the common system as- 
sociated with the central rapidity, only the second and fourth 
quadrants in Fig. [7] would be filled. 




A FC [deg] 



FIG. 7. (Color online) The 2-dimensional distribution plot 
of the relative torque angles Afc and Abc, for centrality 
50 — 60%, space-time rapidity 7711 = 2.5. The corresponding 
correlation coefficient is pfcb = —0.61. 
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FIG. 8. Covariance of the Afc and Abc angles plotted as a 
function of the space-time rapidity 7711 for various centralities. 



III. HYDRODYNAMICS 

The intermediate evolution of the dense system formed 
in relativistic heavy-ion collisions is believed to be gov- 
erned by hydrodynamics [15l . |23| . |24j . The hydrodynamic 
expansion of the fireball generates a collective velocity 
field of the hot fluid. The direction of the acceleration of 
the fluid element is given by pressure gradients. This way 
the azimuthal eccentricity of the fireball is transformed 
into the elliptic flow [25], the triangular deformation into 
the triangular flow [221 ]. and the source tilt into the di- 
rected flow of the final hadrons Q . A similar mechanism 
generates, on event-by-event basis, a torqued transverse 
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FIG. 9. Same as in Fig. [8] for the rank-3 angles (triangular 
flow). 



velocity field at different rapidities. 

The discussion of the torque of the fireball in the pre- 
ceding sections concerned the earliest stage of the col- 
lision, described within the wounded-nucleon approach. 
That stage, essentially, prepares the initial condition for 
the subsequent phases of the evolution. It is necessary 
to check that the signatures of the torque fluctuations 
survive these later stages, such that measurable effects 
can be detected in experiment. Here we investigate the 
behavior of the rank-2 torque under the hydrodynamic 
evolution, as the rank-3 case is expected to behave sim- 
ilarly. A rotation of the density in the transverse plane 
in the torqued fireball scenario, would generate a rotated 
fluid velocity field. This rotation of the transverse ve- 
locity field at each space-time rapidity would lead to a 
rotation of the reconstructed reaction plane for particles 
emitted in the corresponding rapidity range. 

We apply the 3+1 hydrodynamic evolution of the per- 
fect fluid with a realistic equation of state, implemented 
for the first time in [26j . This approach is capable of uni- 
formly describing the main experimental features, such 
as the transverse momentum spectra, v-z. as well as the 
Hanbury Brown-Twiss correlation radii [24j . Extension 
to 3+1 dimensions allows to describe also the spectra at 
non-central rapidities and the directed flow v\ [1, [ll| . 

To demonstrate the effects of the torqued fireball on 
final particle spectra, we generate the hydrodynamic evo- 
lution for the sample value of the impact parameter, 
b = 6.6 fm, which corresponds to the centrality c = 20- 
25% [I?} . The initial energy density e is taken in the form 
given in Eq. @, with the wounded nucleon and and bi- 
nary collision densities calculated in the Glauber model. 
The details and the parameters of the distribution can be 
found in 0, [ll| . The initial proper time for the evolution 
is To = 0.25 fm/c. To study the torque effect, instead of 
doing tedious event- by- event hydrodynamic simulations 
with a whole distribution of torque angles as in Fig. |H 
we perform one simulation where the densities of the for- 
ward and backward going wounded nucleons are rotated 
in opposite directions by a fixed angle of 5° which is a 
value corresponding to the rms width of the distribution 
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FIG. 10. (Color online) The dependence of the torque an- 
gle of the fluid velocity field on space-time rapidity after the 
3+1-dimensional hydrodynamics of Ref. (solid, dotted and 
dashed lines). Subsequent curves are for different evolution 
times. 



of Afb- Thus the initial energy density ©, which is the 
starting point of our hydrodynamics, becomes 

^(V\\,x,y) = (1 - a)[p + {Rx,Ry)f + (r)\\) 

+p- (R T x, R T y)f^ (77|, )] + ap hin (x, y)f(v\\)- (11) 

The operator R rotates the density of forward going 
wounded nucleons by the fixed angle, while operator R T 
rotates the backward going wounded nucleons in the op- 
posite direction. The density of binary collisions is not 
rotated, as its emission component is symmetric in 7711 . 

During the hydrodynamic evolution the direction of 
the transverse flow is determined by the orientation of 
the fireball density at a given space-time rapidity. The 
torque angle of the fluid, ^(ffli), is determined by the 
requirement that in the frame defined by the principal 
axes of the transverse flow we have (T xy }(r)n) = 0, where 
(!),„) (?7[|) = J dxdyT^(r]\^x,y) are the components of 
the energy-momentum tensor averaged over the trans- 
verse plane. In Fig. [TU] we show the evolution of the 
torque angle of the fluid velocity field in the hydrody- 
namic calculation. We present the angle as a function of 
the space-time rapidity after a hydrodynamic evolution 
lasting 1, 3 or 6 fm/c. Due to the longitudinal push, 
the effect decreases somewhat as the evolution time in- 
creases, but this quenching is very small. Therefore the 
torque effect survives the hydrodynamic phase with an 
almost unchanged magnitude. A sizable twist of the col- 
lective velocity field at the freeze-out should give a twist 
of the distribution of the emitted particles. 

We note that the space-time rapidity ?y|| is not equal to 
the fluid rapidity, as in a boost-invariant model. More- 
over, the transverse momenta and rapidities of final 
hadrons include a thermal component besides the col- 
lective velocity. These effects are addressed in the next 
section, through the use of a realistic model of statistical 
hadronization. 
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FIG. 11. (Color online) The event-by-event distribution ol k6 
for Vk = 5% for several values of the event multiplicity n: 600 
(solid), 100 (dashed), and 20 (dotted). 



IV. HADRONIZATION AND EXPERIMENTAL 
MEASURES 

The final question, of key practical importance for the 
whole idea, is how to observe the torque of the fire- 
ball from the data consisting of momenta of detected 
particles. In our model approach we adopt the statis- 
tical hadronization picture (28[, where hadrons (stable 
and resonances 29], which subsequently decay) are pro- 
duced at freeze-out according to the Frye-Cooper for- 
malism [lO], with a freeze-out temperature of 150 MeV 
Hill. The difficulty lies in the fact that the finiteness 
of the number of particles produced in each event causes 
sizable fluctuations of the principal axes (or the event 
plane) as determined from the transverse momenta. We 
denote the relative angle between the event-plane axis 
and the fireball spatial principal axis \I/ as 9. Despite 
this difficulty, as we shall see, one can propose measures 
of even-by-event torque fluctuations that should be pos- 
sible to be observed in the RHIC data. 

Throughout this section r\ denotes the momentum 
pseudorapidity of produced particles, in distinction of the 
space-time rapidity m of the preceding parts. 



Fluctuations of principal axes from statistical 
hadronization 



To appreciate the phenomenon of the fluctuation of 9, 
let us recall formulas from [l(| concerning the fluctua- 
tion of the eccentricity and the principal axes due to the 
finite number of (independent) particles. The angle 9 
(dependent on the rank k) is defined in each event as 



tan(fc0) 



ELi sin ( fc ^) 
E"=i cos(kfa) ' 



(12) 



where all angles are measured with respect to the angle 
of the principal axis of the fireball, &( k \ at a selected r\ 



window. From Appendix C of Ref. [Ioj |. which derives the 
results based on the central limit theorem, we find that 
(at a fixed multiplicity n) the event-by-event distribution 
of 9 is given by 



f(k9) 



2ns 3 
sgn[cos(fc0)]erf 



{ V^v k cos(k9)e nv * ™s 2 (fc0)/( s V) 
'nufcl cos(k9)\ ' 



sw 



sw 



(13) 



where the short-hand notation is 



w = \ll-2v 2 k , 



s=Jl- 2v 2 k sui\k9). 



(14) 



This distribution is a function of the product k9 and 
Vky/nj therefore we expect universality with respect to 
the rank fc, as well as scaling with Vk^/n- Both larger 
multiplicity and larger Vk reduce the fluctuations of the 
angle 9. 

In Fig.[TT]we plot the distribution J(k9) for the typical 
value of the flow coefficient for the elliptic flow (k = 2) 
with v k = 5% for three values of the multiplicity n. We 
note a large width of the distributions, of the order of 30° 
for k = 2, or even more as n or vi decrease. This spread 
will have the tendency of washing out the smaller torque 
angle, unless a suitable method, sensitive to differences 
of angles in the same pseudorapidity bin, is used. 



B. Cumulants 

To achieve the goal of observing the torque, similarly 
to analyses of flow [13, HU, we consider cumulants. In 
the simplest case of the two-particle cumulant we may 
take 



y —yy> 

Advents ^ n F n B f-f f-f 

events %— 1 j — 1 



} ik(4>i-<j>j) 



(15) 



with k denoting the Fourier rank and fa (fa) being the az- 
imuthal angles of particles emitted in the forward (back- 
ward) rj windows in a selected centrality class. The quan- 
tities np and riB are the corresponding multiplicities. 
The measure is averaged over events, with the number 
of events in the sample equal to iV cvcnts . 

When no correlations between particles are present, 
the distribution function of n particles is the product of 
one-body distributions. The one-body distribution can 
then be written in the form 



f(fa =v + 2J2 Vk cos[A(0 -¥<*>)]. 



k=l 



Then 



b)\ = 



V k ,FVk,B COs(fcA FB )) e 



(16) 



(17) 
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where (-) cvcnts indicates the averaging over events. Non- 
flow contributions modify the right-hand side at the level 
1/n, where n is the effective multiplicity of particles in 
the event. These effects, hard to estimate, include res- 
onance decays, conservation laws, Bose-Einstein correla- 
tions of identical particles, short-range correlations, etc. 
The influence of resonance decays will be analyzed via 
simulations below. 

Since we are interested in measuring the average 
cos[n(\&F — ^b)]) we need to divide Eq. (|TT|) by Vk,FVk,B- 
We can do it, for instance, by evaluating the ratio of cu- 
mulants, defined as 

/ e ik(<p F -4> B )\ 

cos(kA FB ) {2} = _=L 1 = = 

(cos(fcA FB )) event8 + norrflow. (18) 

One may also use higher-order cumulants to generate 
measures of the torque. For example, the ratio of four- 
particle cumulants, yields 

l e ik[(4>F,1+<t>F,2) — {<t>B,l+<t>B,2)]\ 

cos(2fcA f B ) {4} = , — r - 7 - ; — Jr = 

(cos(2fcA FB )) evelrts + nonflow (19) 

The practical issue in this kind of studies is the influ- 
ence of the non-flow contributions on the result. 



C. THERMINATOR simulations 

The precise estimate of the non-flow effects in the for- 
mulas of the previous section is not easy. To obtain a 
realistic estimate of the influence of resonance decays, 
we have run THERMINATOR [Hj], generating 100000 events 
(in one centrality class) from a fireball with a torqued 
hypersurface, resulting from running the 3+1 perfect hy- 
drodynamics Q on the torqued initial condition, as de- 
scribed in Sect. IIIII In Figs. [T2]we present the results for 
centrality c = 20 — 25%, where, at freeze-out, the torque 
angles in the forward and backward directions differ by 
~ 8°, as in Fig. [TO] 

First, we check if THERMINATOR is capable of reproduc- 
ing the (fixed) input torque angle. For this purpose we 
take into account primordial particles (those created at 
the freeze-out hypersurface) only, thus disregarding reso- 
nance decays. The result, shown in Fig. [T2]with squares, 
shows a nice agreement between the cumulant measures 
(I18I19P and the functions cos(2Afb) and cos(AAfb), 
shown with lines, evaluated directly from the fireball 
torque angle shown in Fig. [TU] (we take here the case of 
the evolution time equal of 6 fm/c). For comparison, the 
triangles indicate the result of the calculation without 
the torque, Afb — 0. The agreement shows that with- 
out the non-flow contribution from resonance decays the 
accumulated statistics of 100000 THERMINATOR events per 
centrality class is sufficient to see the torque effect. 



The error bars in figures of this section correspond 
to the statistical errors of our Monte Carlo simulation. 
These errors increase fast with rj, as the number of pro- 
duced particles decreases rapidly above \r)\ ~ 3. In exper- 
imental data samples, which at RHIC have a very large 
statistics, these errors should be significantly smaller. 

Next, we show the realistic case, accounting for all 
charged pions, kaons, protons, and antiprotons, includ- 
ing those coming from resonance decays, in the deter- 
mination of the principal axes. We set the transverse 
momentum in the window 0.45 GeV < pt < 3 GeV. 
This limits somewhat the contribution of resonance de- 
cays, which populate predominantly the softer part of 
the momentum spectrum. Also, higher px leads to a 
larger elliptic flow coefficient, which reduces the statisti- 
cal noise (cf. Eq. (fT3"|) and its discussion), thus is advan- 
tageous for our calculation with a limited-size sample. 
The result of our simulation is shown in Fig. 1131 with the 
squares corresponding to the calculation with the torque 
as in Fig. [TU1 and the triangles giving the base-line re- 
sults with no torque. We note that both cos(2A F s){2} 
and cos(4Afb){4} pick up the non-flow effects, as the 
triangles are shifted down from unity. 

We see in Fig. [T3] that the non-flow effect from res- 
onance decays depends weakly on ry. This is to be ex- 
pected, as the effect comes from correlations of the two 
particles in the forward pseudorapidity window and the 
two particles in the backward pseudorapidity window. 
Once these windows are sufficiently separated, the effect 
does not depend on the separation rj. We also remark 
that the simulation is much more noisy (not shown in 
figures) for the central events, which is due to the small 
value of V2- Thus the torque effect has the best chance 
of being observed for the mid-central or mid-peripheral 
centrality classes. Thus, if a decrease of the proposed cu- 
mulant measures with 77 is observed in real data, it would 
hint to the torque effect. 

In a complete simulation, an average should be taken of 
the observables over the torque angle distribution. With 
our calculation, taking the torque angle corresponding to 
the rms width of the angle distribution, we get a correct 
estimate for the averages of quantities of the order of the 
square of Afb (Eqs. IT5I IT"9"|) . which is sufficient. 



D. Forward-backward-central correlations 

On may also use correlation measures based on three 
rapidity windows, forward, backward, and central. One 
possibility is the following combination: 

AfbcH} = 

/ e i2[(<t>F-<t>C,l)-(<t>B-</>0,2)]\ _ ( e i2[(4>F-4>C,l) + (<l>B-4>C,2)]-j 

V2,FV2,BV2,C 

= (2 sin(2A FC ) sin(2A BC )) cvonts + nonflow. (20) 
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FIG. 12. (Color online) Cumulant measures of the torque 
obtained for the Au+Au collisions at c = 20 — 25% with the 
primordial particles only (i.e. with no resonance decays), plot- 
ted as functions of pseudorapidity. Triangles correspond to no 
torque, squares to the torque of Fig. [TO] The solid line repre- 
sents evaluation directly from the fireball torque angle shown 
in Fig. 1101 The agreement of the line and the squares shows 
that the statistics is sufficient to detect the torque effect. The 
r\ windows have the width of one unit. The error bars corre- 
spond to the statistical errors of the THERMINAT0R simulation. 
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FIG. 13. (Color online) Same as Fig. [12] obtained from 
all charged pions, kaons, protons, and antiprotons, with 
450 MeV < pr < 3 GeV. The departure of the triangles, cor- 
responding to no torque, from unity displays the non-flow con- 
tribution due to resonance decays. The torqued case (squares) 
is shifted from the case without the torque. 



V. CONCLUSIONS 



In conclusion, we summarize our results: 



For small torque angles this measure reduces to the co- 
variance ([9]), namely 

AfbcI^} ~ 8cov(A FC , A BC ) + nonflow. (21) 

The results of the THERMINAT0R simulations for A FB c 
are shown in Fig. 1141 The conclusions are similar as for 
the previously discussed measures. With primordial par- 
ticles (no non-flow effects, Fig. Q3Ja)) we reproduce the 
expected behavior. The squares are very close to the 
solid line, showing 8cov F bc (cf. Fig- EJ) - F° r the realistic 
case of all charged pions, kaons, protons, and antiprotons 
(Fig. flir b)). we observe the non-flow effect from the res- 
onance decays. We note that the cases of the torqued 
and untorqued fireballs are qualitatively different, with 
the former giving positive, and the latter negative values 
of A FB c- 



1. The space-time rapidity emission profile, where the 
wounded nucleons emit predominantly in the direc- 
tion of their motion, combined with the statistical 
fluctuations of the source densities in the transverse 
plane, lead to event-by event torqued fireballs. 

2. The standard deviation of the relative torque an- 
gle between the forward (7711 ~ 3) and backward 
(77H ~ —3) regions varies from 20° for the most cen- 
tral collisions to 10° for the mid-central and mid- 
peripheral Au+Au collisions at the highest RHIC 
energies. 

3. This initial torque is transformed, via hydrodynam- 
ics, into the torque of the transverse collective flow 
of the fluid, and subsequently into the torque of the 
principal axes of the transverse-momentum distri- 
butions of the detected particles. 

4. Statistical measures based on cumulants containing 
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FIG. 14. The measure Apse obtained from for the primor- 
dial particles (a) and all charged pions, kaons, protons, and 
antiprotons (b). Triangles correspond to no torque, squares 
to the torque of Fig. 1101 The solid line shows 8cov fbc of 
Eq. (O. The r\ windows have the width of one unit. The error 
bars correspond to the statistical errors of the THERMINATOR 
simulation. 



particles in different pseudorapidity bins should be 
useful in detecting the torque effect experimentally. 

5. The non-flow corrections can be sizable, but should 



not overshadow the effect. In particular this is 
the case of the effects of resonance decays, esti- 
mated realistically with THERMINATOR Monte Carlo 
simulations. We find that there is a clear differ- 
ence between the behavior of the proposed cumu- 
lant measures for the untorqued and torqued cases. 
Therefore the torque fluctuations should possibly 
be observed in the high-statistics RHIC data by 
the PHOBOS and STAR collaborations. 

6. Since the statistical noise increases as the product 
of the particle multiplicity and the Vk flow coef- 
ficient, the best chance for looking for the torque 
effect is in the mid-central or mid-peripheral cen- 
trality classes, such as c = 20 — 30%, and with the 
exclusion of the soft-momentum hadrons. 

7. Based on our calculations, the torque is expected 
to have a similar size for the elliptic flow and tri- 
angular flow. 

8. Other effects that may influence the torque ef- 
fect, such as the conservation of momentum [32J, 
angular- momentum [33[ , or charge [!, H3| , as well 
as other sources of correlations and fluctuations, 
should be incorporated in future studies. 

9. Finally, the torque fluctuations in other two-source 
or multi-source models 35] should be investigated. 
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